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VARIATION  IN  THE  AMPLITUDE  OF  PERTURBATIONS  ON  THE  INNER 
SURFACE  OF  AN  IMPLODING  SHELL  DURING  THE  COASTING  PHASE 

Introduction 

In  most  inertial  fusion  target  designs  the  pellet  shell  goes  through  three  stages  during  the 
implosion:  acceleration,  coasting  and  deceleration.  First  the  shell  is  accelerated  inward, 
with  Rayleigh-Taylor  unstable  growth  at  the  outer  ablation  surface  and  at  any  layer  in¬ 
terface  where  a  lighter  fluid  is  pushing  on  a  heavier  material.  Then  after  the  shell  reaches 
nearly  its  final  velocity,  it  coasts  inward  with  little  or  no  acceleration.  With  slight  over¬ 
simplification  we  can  view  this  phase  as  ballistic.  And  finally,  the  shell  decelerates  as  it 
compresses  the  central  gas  ignitor.  This  third  phase  is  also  Rayleigh-Taylor  unstable,  on 
the  inner  surface  of  the  shell.  For  a  large  high-gain  pellet,  the  first  phase  reduces  the  pellet 
radius  from  about  3000  microns  to  about  1500  microns.  The  second  phase  reduces  the 
radius  further  to  about  150  microns.  As  the  radius  decreases  tenfold,  the  mass  density 
and  the  shell  thickness  each  increase  about  tenfold.  The  third  phase  reduces  the  radius 
further  to  about  75  microns. 

During  the  first  (acceleration)  phase,  Rayleigh-Taylor  perturbations  on  the  outer  sur¬ 
face  drop  off  radially  roughly  as  exp (—lx/R),  where  l  is  the  mode  number  of  the  perturba¬ 
tion,  R  is  the  shell  radius,  and  x  measures  the  distance  inward  from  the  unstable  surface. 
For  many  mode  numbers  of  interest,  the  shell  is  initially  so  thin  that  the  inner  and  outer 
surfaces  have  about  the  same  perturbation  amplitude.  These  inner-surface  pertubations 
can  then  serve  as  a  finite-amplitude  source  for  the  Rayleigh-Taylor  mode  during  the  decel¬ 
eration  phase.  The  purpose  of  the  present  note  is  to  study  the  evolution  of  the  inner-surface 
perturbations  during  the  ballistic  phase  and  to  examine  the  role  of  convergence. 

One  possibility  during  the  ballistic  phase  is  that  the  perturbations  evolve  so  that  the 
mass  of  any  protuberance  or  bulge  is  conserved.  After  all,  the  Rayleigh-Taylor  phenomenon 
moves  mass  from  valleys  into  bulges.  If  the  Rayleigh-Taylor  growth  is  negligible,  one  might 
guess  that  the  mass  in  a  radial  sector  is  conserved  according  to  pR2$  —  const,  where  £  is 
the  amplitude  of  an  inner-surface  perturbation.  We  find  that  the  perturbations  actually 
evolve  quite  differently.  As  the  shell  thickens  and  the  inner  and  outer  surfaces  decouple, 
“waves”  can  propagate  on  the  inner  surface  of  the  shell.  We  will  show  that  the  action 
associated  with  these  waves  is  an  adiabatic  invariant. 

A  realistic  analytic  model  would  use  compressible  surface  perturbations  on  a  com¬ 
pressible  imploding  shell.  Hore,  in  order  to  obtain  analytically  tractable  equations,  we 
assume  that  the  perturbations  are  irrotational  and  incompressible.  This  is  a  plausible 
assumption,  based  on  other  studies  of  Rayleigh-Taylor  instability  in  compressible  and  in¬ 
compressible  media.  For  simplicity  we  will  also  assume  that  the  zeroth-order  fluid  motion 
is  incompressible.  This  second  assumption  is  not  a  good  one  in  inertial  confinement  fusion, 


but  with  this  model  we  can  rigorously  derive  our  adiabatic  invariant  from  the  fluid  equa¬ 
tions,  study  its  properties,  and  gain  some  physical  insight.  We  suggest  that  the  adiabatic 
invariant  may  be  useful  in  the  coasting  phase  even  for  compressible  media. 

In  order  to  motivate  this  approach  it  is  useful  to  illustrate  the  physical  ideas  by  means 
of  the  example  of  a  pendulum  whose  length  is  slowly  changing  (the  “Rayleigh  pendulum”). 
The  total  energy  of  a  simple  pendulum  describing  small  oscillations  of  amplitude  9  in  a 
uniform  gravitational  field  with  acceleration  g  can  be  written  as 

H  =  \-mL2{92  +  g62  /  L)  =  ]-(p$2 /mL2  +  mgLO 2)  =  \-mgLQ2,  (1) 

z  z  z 

where  L  is  the  length,  m  is  the  mass,  the  dot  indicates  a  time  derivative,  and  we  have 
introduced  p$  =  mL29 ,  the  momentum  canonically  conjugate  to  9.  When  H  is  constant 
the  motion  of  the  pendulum  is  given  by  9  =  Qsinu;(t  —  t o),  where  0  is  the  amplitude  and 
u  =  ( g/L j1/2  is  the  frequency  of  the  oscillations. 

If  now  the  gravitational  acceleration  g,  the  mass  m,  or  the  length  L  is  changed  slowly, 
i.e.,  over  a  time  At  such  that  uAt  »  1,  then  H  is  no  longer  conserved  but  the  motion  is 
still  approximately  harmonic.  By  a  well-known  theorem,  we  can  invoke  the  constancy  of 
the  action  §  pg  d9  =  irmu>L2Q2 ,  where  the  integral  is  carried  out  over  a  complete  cycle  of 
the  motion,  to  compute  the  change  in  0.  For  example,  if  it  is  the  length  which  is  varied,  0 
changes  according  to  I,3/2©2  =  const.  This  result  is  of  course  identical  with  that  obtained 
by  direct  calculation  of  the  work  done  on  the  pendulum  by  the  time-averaged  gravitational 
and  centrifugal  forces. 

On  the  basis  of  this  example  we  anticipate  that  the  action  of  a  stable  small-amplitude 
oscillation  on  the  surface  of  an  imploding  spherical  shell  will  be  a  good  adiabatic  invariant 
if  the  implosion  time  is  long  compared  with  the  period  of  the  oscillation.  From  dimensional 
considerations,  we  expect  the  action  to  be  proportional  to  the  product  of  the  square  of  the 
perturbation  amplitude  the  mode  frequency  u>,  and  an  effective  mass.  Because  the  inner 
surface  of  a  coasting  shell  accelerates  inward  on  account  of  convergence  effects,  waves  on 
the  inner  surface  are  stable  and  satisfy  the  dispersion  relation  u)  «  (—IRi/Ri)1/2,  where  l  is 
the  mode  number.  If  the  eigenfunctions  satisfy  Poisson’s  equation,  they  will  decay  toward 
the  interior  of  the  shell  over  a  characteristic  distance  of  the  order  of  the  inverse  wavenumber 
R\jl.  The  effective  mass  will  thus  be  approximately  {AnR\2){R\jl)p  —  \tt R\Z  pjl,  and  the 
adiabatic  invariant  will  be 

4.nRi3p$2u/l  «  const.  (2) 

Consequently  we  expect  that  for  a  given  l  the  amplitude  of  a  surface  perturbation  will  vary 
approximately  as  (— R\Z  R\)~1^ . 

We  begin  our  investigation  by  analyzing  the  motion  of  a  coasting  spherical  shell.  To 
determine  in  detail  how  a  surface  wave  evolves,  we  then  linearize  the  ideal  equations  of 
motion  for  a  shell  of  uniform  density  and  solve  them  under  the  assumption  that  the  flow 
is  incompressible  and  irrotational.  In  the  following  two  sections  we  look  at  the  thin-shell 
and  thick-shell  limits.  When  the  shell  is  thin  the  perturbations  grow  linearly  in  time.  For 
thick  shells  the  perturbations  localized  on  the  inner  surface  decouple  from  those  on  the 
outer.  We  show  that  for  certain  forms  of  the  motion  the  linearized  equations  can  be  solved 
exactly.  In  the  limit  of  thick  shells  and  large  mode  number  l  we  carry  out  a  WKB  analysis 


of  the  equation  for  perturbations  on  the  inner  surface  and  recover  the  adiabatic  invariant 
(2).  To  check  our  analysis  we  then  solve  the  exact  linearized  equations  numerically.  In  the 
last  section  we  discuss  these  results  and  offer  our  conclusions. 


Coasting  Shell 

The  total  mass  and  total  energy  of  the  imploding  shell  are  conserved.  The  mass  is 

M  =  4tt  [**  p(R)R 2  dR  =  ^p{R23  -Ri3)  =  ^ pRo 3,  (3) 

J  R1  *5  O 

where  R\  and  R2  are  the  inner  and  outer  radii  of  the  shell  at  time  £,  p(R)  is  the  density 
as  a  function  of  radius,  p  is  the  average  density,  and  Rq  is  the  radius  of  the  fully  collapsed 
shell.  Because  the  shell  is  assumed  to  be  incompressible,  the  volume  between  Ri  and  an 
intermediate  fluid  element  moving  with  time-dependent  radius  R  is  conserved  and  we  can 
differentiate  to  get  the  velocity  V(iE,t): 


R2V  =  R2R  =  R^Rl 


Hence  the  total  energy  is 


r  Ri 

W  =  2n  p(R)R2R2  dR 
**2npR12R14  J*3  = 


R 1  jp 

—  =  2irpRi2Ri3(l  —  Ri/R2). 


For  a  uniform  shell,  p{R)  =  p  and  Eq.  (5)  is  exact.  We  assume  this  in  what  follows. 

The  characteristic  time  associated  with  the  implosion  is  {MRq2 /2W)1/2 .  It  is  conve¬ 
nient  to  choose  the  zero  of  time  to  correspond  to  R\  =  0.  By  differentiating  Eq.  (5)  and 
using  Eq.  (4)  with  R  =  R2  we  find  the  equation  of  motion  of  the  unperturbed  shell, 


Ri  = 


AirpRi4 


o  R\  (Rxy 

3+2ir  +  hr  ; 


R2  satisfies  a  similar  equation.  Both  are  easily  solved  numerically  (Fig.  1).  Note  that 
Ri  <  0  <  R2,  so  that  perturbations  on  both  surfaces  are  stable.  The  quantity  in  square 
brackets  in  Eq.  (6)  changes  by  only  a  factor  of  two  as  Ri/R2  varies  between  0  and  1,  so 
to  a  good  approximation  we  have  cx  R\~4,  whence  Ri  oc  (— t)2/5  and 


Ri/R\  oc  R\  5  oc  (— t)' 


3 


Linearized  Equations 

For  an  incompressible  fluid  the  linearized  continuity  equation  can  be  written 


V  •  v  =  0,  (8) 

where  v  is  the  first-order  velocity.  (Here  and  throughout  this  paper,  upper-case  variables 
are  associated  with  the  unperturbed  motion  and  lower-case  variables  with  the  perturbed 
motion.)  If  the  perturbations  are  assumed  to  be  irrotational,  we  can  introduce  a  scalar 
potential  by 

v  =  V<f>.  (9) 

This  potential  satisfies  the  Laplace  equation, 

V2<f>  =  0.  (10) 

The  solution  of  Eq.  (10)  can  be  written 


<£(R,0  =  Y,[Riuim(t)(Ri/R)l+1  +  R2Ulim{t)(R/R2y}Yltrn(0,ip),  (11) 

l,m 

where  Uj,  j  =  1,  2,  plays  the  role  of  a  time-dependent  amplitude  of  the  velocity  pertur¬ 
bation  associated  with  R  =  Rj ;  Y;)m  is  the  spherical  harmonic  of  order  Z,m;  and  the  sum 
is  taken  over  all  meridional  and  zonal  mode  numbers  such  that  0  <  |m|  <  Z  <  oo.  Since 
the  perturbed  equations  are  linear,  in  what  follows  we  specialize  to  a  particular  mode  and 
omit  the  superscripts  Z,m  whenever  doing  so  will  cause  no  confusion. 

The  velocity  associated  with  this  mode  is  found  by  taking  the  gradient  of  Eq.  (ll): 

v(R,<)  =  {tt1(J21/E),+2[-er(/  +  l)Yt,m  +  RVYitm] 

+  u2(R/R2)'-1  [erlYi,m  +  EVy/>m]  },  U2) 

where  er  is  the  unit  vector  in  the  radial  outward  direction.  The  linearized  momentum 
equation  can  be  written 

+  Vp  =  0.  (13) 

To  find  the  first-order  pressure  p(R,Z)  we  employ  the  Bernoulli  form  of  Eq.  (13),  obtained 
by  combining  Eq.  (9)  with  Eq.  (13)  and  integrating: 

P  +  Vv^j  +P  =  /(*)•  (H) 


dv 

dt 


+  V(V-v) 


The  function  of  t  on  the  right  side  can  be  made  to  vanish  by  proper  choice  of  ($>.  Substituting 
Eqs.  (ll)  and  (12)  in  Eq.  (14),  we  obtain 


P 

P 


+ 


R\U\  f  (Z  -)-  2) R\U\  —  (Z  +  l)/2] uj  (  — 


}h 

R 


i ;  i 


R2U2  —  (l  —  l)/Z2u2  -f  IR2U2 


R2 

R 


R 

Ri 


Yi 


l ,m  ■ 


(15) 


The  boundary  conditions  needed  to  complete  the  specification  of  the  problem  are  the 
kinematic 


and  dynamic 


»*«*(&),. 

p,+(^),6=° 


conditions  for  j  =  1,  2,  where  ft  is  the  amplitude  of  the  first-order  displacement  of  the  jth 
surface.  Equations  (16)  and  (17)  are  valid  regardless  of  the  dynamics  of  the  basic  state; 
i.e.,  they  hold  for  any  form  of  Ri(t). 

Using  Eq.  (4)  and  the  unperturbed  equation  of  motion  in  the  form 


pR  =  -  VP, 


we  can  rewrite  Eqs.  (16)  and  (17)  as 


vrj  —  ft  +  2  Rj  ft  j  Rj 


Pi  —  pRjCj  >  (20) 

j  =  1,  2.  Setting  R  =  Rj  in  the  radial  component  of  Eq.  (12)  and  substituting  the  result 
in  (19)  yields  two  linear  ordinary  differential  equations  relating  ft  and  ft  to  «i  and  u2, 
which  can  be  written  as 

ft  +  2(Ri/Ri)ft  =  — (/ +  l)ui  + lU2  (21) 


ft  +  2(i?2/i?2)ft  =  —{l  +  l)Al+2u\  +  Zu2, 
where  A  =  Ri/R2.  Combining  Eq.  (15)  with  Eq.  (20),  we  have 

(RjUi)  +  ji?2u2  +  \lA  3  —  l  +  l]  P2U2  j  A1  =  —  .Rift 


aria 

+  [/  +  2  —  (Z  +  l)A3]  R\U\  |  Al+1  +  (R2U2)  =  —R-2^2-  (24) 

Equations  (2l)-(24)  constitute  a  fourth-order  system  of  linear  ordinary  differential 
equations  in  ft,  ft,  ui,  and  u2.  Before  solving  them  numerically  we  analyze  them  in  the 
limits  in  which  the  shell  is  thin  or  thick  compared  with  the  radius. 


Thin-Shell  Limit 

When  the  shell  thickness  r  =  i?2  —  Ri  is  much  less  than  either  R\  or  R2,  we  have  to  a 
good  approximation 

r  Rq 

R  =  3R*'  ^ 

where  R  =  +  R2).  Hence  we  can  expand  the  perturbation  amplitudes 

and  the  coefficients  in  the  velocity  potential  Uj  =  u®  +  uj,  j  =  1,  2,  where 

Cy/Cy  ~  u}/uj  ~  r/R  <  1.  (26) 

By  Eq.  (6)  and  its  analog  for  I?  2  the  accelerations  Rj  are  first-order  quantitities,  so  to 
lowest  order  Eqs.  (23)  and  (24)  both  yield 


[i?(u°  +  u°)]'=  0, 


whence  we  obtain 


u°2  =  - ul  +  C/R , 

where  C  is  constant.  Hence  Eqs.  (21)  and  (22)  in  lowest  order  become 

+  2  — =  ~i21  +  !)ui  +  -J£ 

and 

;0  ,  n^.O  _  1>,r°  (^  +  l)C 

£2  +  2R _  +  1)“2 - ^ - • 

In  first  order  Eqs.  (23)  and  (24)  yield 

[R{u\  +  «$)]  -  \irui)  +  ^  [(^2)  -  3J2tt§  =  23^-4 

and 

+  “2)]  -  ^(rui)  +  ^2)  - 

Subtracting  Eq.  (32)  from  Eq.  (31),  we  obtain 


,  ,  .  ,0  R  o  SH'fP  3lCRr 

(21  +  1)  rti?  —  2— ru?  = - - — . 

v  '  1  I?  1  7rpI24  ■npR2 


Differentiating  Eq.  (25)  twice  with  respect  to  t,  we  find  that  r  =  —2 (R/R)r  and 
r  —  6 (R/R)2r.  If  we  multiply  Eq.  (29)  by  r,  differentiate,  substitute  Eq.  (33)  in  the 
right-hand  side,  and  use  the  expressions  for  r  and  its  derivatives,  we  find  that  satisfies 


R2  npR03R2 


??  =  0, 


where  we  have  used  the  thin-shell  limit  of  Eq.  (5).  Going  through  similar  operations  using 
Eq.  (30)  instead  of  Eq.  (29),  we  find  that  the  second  derivative  of  $$  likewise  vanishes.  Thus 
both  £1  and  f 2  vary  linearly  with  time.  This  simple  result  should  come  as  no  surprise:  when 
the  unperturbed  motion  of  the  shell  is  uniform,  so  is  that  of  the  perturbations. 


Thick-Shell  Limit 

Equations  (2l)-(24)  also  simplify  considerably  in  the  thick-shell  limit.  If  A1  <C  1,  the 
perturbations  localized  on  the  inner  and  outer  surfaces  decouple  and  satisfy 


-lk«” 


respectively.  We  are  particularly  interested  in  the  behavior  of  the  former.  Here  we  discuss 
several  examples  for  which  Eq.  (35)  can  be  solved  exactly. 

(t)  Coasting  shell.  If  we  substitute  Eq.  (7)  in  Eq.  (35),  we  find  that  a  (— t)a ,  where 

a  =  ±— n/25^2 41.  (37) 
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For  /  >  1,  a  is  complex  with  real  part  equal  to  —1/10.  Hence  as  t  — ♦  0,  we  have  |fj|  ~ 
(— f)-1/10  ~  i?!-1/4  by  Eq.  (7).  This  divergence  is  much  weaker  than  that  (~  i?i-2) 
predicted  by  the  “constant-mass  perturbation”  model. 

(»)  Constant  acceleration.  The  linearized  equations  derived  above  hold  for  any  Ri(t),  even 
one  which  does  not  conserve  the  energy  of  the  unperturbed  state.  If  the  shell  undergoes 
constant  acceleration  G  <  0  and  reaches  the  origin  at  t  =  0  with  velocity  U  <  0,  then  its 
radius  as  a  function  of  time  is  given  by 

Rx{t)  =  \Gt2 +  Ut.  (38) 

Evidently  Vx  vanishes  and  R\  attains  its  maximum  at  to  =  —U/G.  Equation  (35)  becomes 

(|Gt2  +  Ut)f i  +  3 {Gt  +  C/)?x  -  (/  -  l)Gft  =  0.  (39) 

If  we  set  x  =  —Gt/2U,  then  z  =  Ci  satisfies  the  hypergeometric  equation1 

x(l  —  x)z"  4-  [c  —  (a  +  6  +  1)1)2'  —  ahz  =  0,  (40) 


0  5:  2  <  2’  w^ere  primes  denote  differentiation  with  respect  to  x  and 

a  —  ^(5  +  \/lf  4"  8/); 

6  =  ^(5  -  v/17  +  8 /); 
c  =  3. 


(41a) 

(416) 

(41c) 


One  solution  of  Eq.  (40)  is  the  hypergeometric  function  zx  =  F(a,6;c;x),  which  equals 
unity  at  x  =  0  and  takes  on  the  value  2v/7T  [r(|a  +  |)T(|6  +  |)]  at  x  =  1/2.  Since  z\ 
is  related  to  the  associated  Legendre  functions  by1 


the  series  for  z\  reduces  to  a  polynomial  in  x  when  (17  +  8/)1/2  is  an  odd  integer,  i.e.,  for 
l  =  1,4,8,13,....  The  second  solution  z-i  is  finite  at  x  =  1/2  but  diverges  at  x  —  0  as 
22  ~  i~2  ~  [—t)~2.  This  suggests  that  the  growth  of  the  inner  surface  perturbations  is 
more  serious  for  an  accelerated  shell  than  for  one  which  coasts. 

(m)  Simple  harmonic  motion.  If  the  shell  radius  varies  with  time  according  to 

Ri(t)  =  Ri  sin(— Clot),  (43) 

where  and  Ho  are  constant,  then  we  can  set  x  =  sin2  flot  and  again  arrive  at  Eq.  (40), 
with 

a  =  j(3  +  \/5  +  4/);  (44a) 

b  =  i(3-  v/5  +  4/);  (44 b) 

c  —  2.  (44c) 

The  two  solutions  behave  similarly  to  those  found  in  case  (tV):  z\  is  finite  for  0  <  x  <  1 
and  reduces  to  a  polynomial  in  x  if  (5  +  4/ ) 1/2  is  an  odd  integer,  i.e.,  for  l  =  1,5,11,19,...,- 
22  diverges  at  t  =  0  as  l/x  ~  (— 1)~2. 

(iv)  Mode  number  1  =  1.  For  l  =  1  Eq.  (35)  becomes 

■Rifi+3.RiCi  =0,  (45) 

which  has  two  types  of  solutions,  satisfying 

=  const  (46) 

and 

o  , 

R\  fi  =  const,  (47) 

respectively.  For  a  thick  coasting  shell,  when  R\  <C  R2,  Eq.  (47)  can  be  rewritten  using 
Eq.  (5)  as 

T7T  ~  const!  (48) 

UJl  1 

which  implies 

•Ri  ~  const.  (49) 

In  general  <a(t)  is  a  superposition  of  these  two  solutions.  The  second  one  blows  up  as 
R\  —*  0,  however,  and  therefore  dominates  for  sufficiently  small  R\.  Consequently  we 
expect  to  diverge  as  R\~ 1//2  in  this  limit.  However,  as  is  well  known,  this  apparent 
divergence  is  merely  a  shift  in  the  location  of  the  origin. 


WKB  Analysis  for  Thick  Shells 

We  look  for  a  rapidly  oscillating  solution  to  Eq.  (35)  in  the  form 


?i  (0  =  z{t)  exp 


(50) 


where  the  argument  of  the  exponential  is  assumed  large  compared  with  unity  (the  WKB 
or  quasiclassical  approximation).  Substituting  in  Eq.  (35),  we  obtain 

— u2  R\z  +  i(3uiRiZ  +  2u)R\z  +  CjR\z )  +  Ri'z  +  3R\Z  +  2Ryz  =  (/  +  \)Riz.  (51) 

In  lowest  order  we  find  an  expression  for  the  square  of  the  mode  frequency, 

u'1(t)  =  -{l  +  l)Ri/Ri.  (52) 

Evidently  the  approximation  is  good  when  /  ^>  1.  In  next  order  we  have 

3uR\z  +  2uR\Z  +  loR\z  =  0,  (53) 

whence 

R\z2oj  =  const.  (54) 

This  is  nothing  other  than  the  adiabatic  invariant  (2).  It  is  interesting  to  note  that  when 
Ri  oc  (  —  t )2/5  (the  form  which  holds  for  a  coasting  shell),  Eq.  (54)  predicts  that  the 
amplitude  of  the  perturbation  varies  as  (i?i)-1/4,  the  exact  answer  found  in  example  (i) 
above. 

The  result  (54)  can  be  reached  by  another,  slightly  different,  approach.  In  the  thick- 
shell  limit,  the  total  kinetic  energy  associated  with  the  part  of  the  perturbation  (11) 
localized  at  the  inner  surface  can  be  represented  as 


Wk 


(55) 


where 


i  f°°  /  ft  \  2i+4 

WKm=\pjR  dRR2(j[)  «i*  "Vl"T 

dft  [-c r(l  +  l)Yt,m  +  RVYi,m]  •  [— er (/  +  l)yi,m  +  RVY,,t 

-IJ-  JJ  dn  [(/  +  i)2  |yJ>m!2  +  R2  |vy,m|2] 


(56) 


hi) 


u 


l,m 


21  +  1 
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■  <- 


-  (- 


v  o  o  * 


ilol  _L  i  ^ 


where  dCl  is  the  differential  of  solid  angle  and  we  have  used  the  fact  that  v  is  real;  the  aster¬ 
isk  (*)  denotes  complex  conjugation.  The  effective  gravitational  energy  of  the  perturbation 
due  to  the  inward  acceleration  of  the  shell  can  be  expressed  as 

WG  =  -pRi2Ri  Jjdn  dqq 

=  — - pRi2Ri  jj  dn  ^!2 

--\pRi2R^\Am]\  jj dn\Yltm\2  (57) 

l,m  *  “ 

I  i,m|2 

=  -2 *pR,*Ri  £  Lj-L 


[In  Eqs.  (56)  and  (57)  we  have  restored  the  superscripts  l,m  labeling  the  coefficients  u,l’m 
and  in  the  summations.]  Combining  these  expressions,  we  see  that  the  total  energy 
associated  with  the  perturbation  can  be  expressed  as  a  sum  over  all  mode  numbers  of  the 
energy  Wl,m  =  W'jf'1  +  WjS71  associated  with  each  (/,m)  mode.  Here 

w>’m  =  2Hf r{(' +  l)K,’T  -  f%l”T} 

2 npRl  f  Id,  2  J,m\  (/  +  l)Rl  I  / ,m  |2  )  /^„1 

-  (/TIj(2/  +  1)  l  ll?dt(Rl  >  -—rT'^  I  I  (58) 

r  (i  +  i)Ri  ,  ‘ 

M  *7l, m9l,m  Ql,mQl,m  » 

where  p  =  27t p/(l  +  1)(2/  +  1) jRj ,  qiiTn  =  -Ri 2 S"i’m ?  and  we  have  used  Eq.  (21)  in  the 
thick-shell  limit  A  —  >  0  to  replace  u^m  with  an  expression  involving  f|’m. 

Equation  (58)  can  be  put  in  Hamiltonian  form.  The  momenta  canonically  conjugate 
to  qitm  are  defined  by  pl  rn  —  pqf  m  for  all  l,m.  With  this  substitution,  Eq.  (58)  becomes 

Wl'm  =  P~1pi,mPlm  +  Ul2qi,mqiymi  (59) 

where  w/2  =  —  (/-f-l)Ri/Rx.  Equation  (35),  Hamilton’s  equation  of  motion  for  qi,m,  follows 
directly  from  Eq.  (59).  The  adiabatic  invariant  (54),  which  is  essentially  / pl  m  dql  m ,  also 
follows  from  Eq.  (59)  by  the  same  theorem  used  earlier  to  calculate  the  adiabatic  invariant 
for  the  simple  pendulum.  (In  fact  the  two  Hamiltonians  arc  formally  almost  identical.) 
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Numerical  Solution 

Equations  (2l)-(24)  can  be  solved  numerically  for  arbitrary  functions  Ri(t),R2(t).  For 
each  choice  of  the  basic  state,  each  value  of  the  mode  number  l  in  general  gives  rise  to 
four  independent  solutions.  An  exhaustive  survey  of  the  behavior  of  surface  perturbations 
is  thus  clearly  impractical.  Instead,  we  limit  ourselves  to  studying  the  evolution  of  a 
representative  sample  of  the  solutions  corresponding  to  a  single  choice  of  basic  state:  the 
coasting  shell  which  served  as  the  original  focus  for  this  investigation. 

We  suppose  that  the  shell  enters  the  coasting  phase  with  small  imperfections,  possibly 
as  a  result  of  instabilities  which  developed  during  the  ablative  drive  stage.  Since  the  shell 
is  very  thin,  these  irregularities  are  likely  to  take  the  form  of  “wrinkles,”  i.e.,  ripples 
which  do  not  alter  the  shell  thickness,  so  we  set  We  can  assume  that  they  are 

initially  quasistatic,  so  that  t*i  =  u2  =  0.  Since  the  perturbed  equations  are  linear,  we  can 
Fourier-analyze  and  set  the  initial  amplitudes  of  each  mode  equal  to  unity  without  loss  of 
generality.  We  follow  the  evolution  of  ten  different  harmonics,  with  mode  numbers  given 
by  1,  2,  5,  10,  20,  50,  100,  200,  500,  and  1000.  The  initial  velocity  (or  energy)  of  the  shell 
determines  the  time  scale  of  both  perturbed  and  unperturbed  equations;  thus  the  only  free 
parameter  in  the  problem  is  the  initial  shell  inner  radius  i?i(£o)-  By  taking  this  to  be  very 
large  compared  with  Ro,  the  radius  of  the  sphere  which  results  when  the  shell  collapses 
completely,  we  can  ensure  that  the  solutions  depend  on  this  choice  only  trivially. 

Our  numerical  results  were  obtained  on  a  VAX  11/780.  The  unperturbed  and  per¬ 
turbed  equations  were  advanced  simultaneously  using  a  fourth-order  Runge-Kutta-Gill 
scheme.  In  our  numerical  solution  we  retained  the  convention  adopted  in  Fig.  1,  according 
to  which  radial  distances  are  measured  in  units  of  R o  and  times  in  units  of  Ro  divided  by 
the  initial  shell  speed.  We  took  the  initial  value  of  the  shell  inner  radius  to  be  equal  to 
100,  corresponding  to  an  initial  time  to  =  —99.28.  This  means  that,  by  Eq.  (25),  the  ratio 
of  the  shell  thickness  to  the  radius  is  r/R  «  3  x  10~7.  This  is  the  amount  by  which  the 
ratio  A  of  the  shell  surface  radii  differs  from  unity.  Since  the  determinant  of  tij  and  u2  in 
Eqs.  (23)-(24)  equals  1  —  A2l+1,  the  equations  are  poorly  conditioned  initially,  and  it  was 
necessary  to  solve  them  in  double  precision  (52-bit  mantissas).  The  time  step  was  chosen 
to  be  St  =  0.01/2i/|i?i|;  with  this  choice  R\  decreased  by  a  factor  of  10s  in  1140  timesteps. 
By  any  realistic  standard,  five  orders  of  magnitude  is  an  enormous  radial  compression. 
The  reason  for  doing  this  was  to  include  in  the  example  both  the  thin-  and  the  thick-shell 
regimes,  and  to  make  sure  that  the  choice  of  initial  radius  did  not  play  an  excessively  large 
role.  The  accuracy  of  the  calculation  was  monitored  by  checking  conservation  of  mass  and 
energy  and  by  comparing  the  numerical  results  with  those  obtained  analytically.  Refining 
St  by  a  factor  of  2  produced  a  negligible  change  in  the  results.  Most  of  the  variation  in  the 
solutions  takes  place  where  R\  deviates  significantly  from  linear  dependence  on  t,  between 
t  —  -1  and  t  =  0  [Fig.  1(6)],  so  we  plotted  all  results  against  log/?!  instead  of  t  or  Rx. 

Figure  2  shows  plots  of  {Ri/R2)1  as  a  function  of  logi?i  for  the  ten  values  of  l  used. 
If  the  transition  from  “thin”  to  “thick”  is  defined  to  occur  when  ( R\/R2 )l  =  0.5,  then  we 
can  approximate  the  transition  radius  by 

R\tkOA1oa.  (60) 

Figure  3  displays  the  linear  growth  (which  appears  curved  because  of  the  logarithmic  scale 
of  the  abscissa)  predicted  by  the  thin-shell  approximation.  In  Fig.  4  the  amplitude  of  the 


perturbation  on  the  inner  surface  obtained  from  the  adiabatic  invariant  (54)  is  plotted, 
also  normalized  to  unity  at  R\  —  100.  Note  that  it  describes  amplification  by  slightly  more 
than  a  factor  of  twenty  over  the  range  of  the  calculation.  As  can  be  seen,  the  thin-shell 
and  the  thick-shell  approximations  cross  near  R\  =  1.25,  i.e.,  in  the  middle  of  the  range 
of  the  transition  radii  defined  by  Eq.  (60). 

Figures  5(a)  and  5(6),  respectively,  show  ft  and  ft  as  functions  of  Rx  for  1  =  1.  Note 
that  the  two  amplitudes  behave  very  differently:  ft  grows  monotonically,  while  ft  reaches 
a  maximum  of  2.9,  then  approaches  unity  (the  initial  value)  as  R\  — *•  0.  The  dashed  curve 
on  Fig.  5(a)  represents  the  approximation  (49)  with  the  constant  set  equal  to  11.2.  As  can 
be  seen,  it  is  reasonably  close  to  the  exact  solution  throughout  the  range  10~3  <  J?i  <  100. 

Figure  6  displays  the  time  dependence  of  the  calculated  values  of  the  perturbation 
amplitudes.  Note  that  the  amplitudes  for  the  inner  (solid  trace)  and  outer  (broken  trace) 
coincide  and  follow  the  thin-shell  approximation  shown  in  Fig.  3  until  Rx  «  R\T.  As 
f?i  ->  0  we  see  that  ft  diverges  while  oscillating  as  a  function  of  log  Ri  with  a  frequency 
that  scales  roughly  as  Z1/2;  in  contrast,  ft  approaches  a  finite  value  of  order  unity.  Perhaps 
the  most  significant  fact  that  emerges  from  these  results  is  that  ft  grows  by  a  factor  of 
about  twenty  while  i?i  is  changing  by  five  orders  of  magnitude. 

To  compare  this  behavior  with  the  prediction  obtained  from  the  conservation  of  action, 
which  is  expected  to  work  best  in  this  limit,  we  have  replotted  ft  in  Fig.  7,  scaled  by  the 
amplitude  shown  in  Fig.  4.  It  is  clear  that  the  adiabatic  invariant  does  an  extremely 
good  job  of  predicting  the  envelope  of  the  oscillations  for  all  l  >  1.  Not  only  is  the  envelope 
of  the  rescaled  amplitude  nearly  constant  through  the  last  three  decades  of  variation  of  Ri, 
but  the  constant  value,  which  is  approximately  equal  to  unity,  is  almost  independent  of  Z. 
Except  for  the  Z  =  1  case  (which  corresponds  merely  to  displacement  of  the  shell  center 
without  deformation),  the  thin-shell  approximation  and  the  adiabatic  invariant  between 
them  thus  describe  essentially  the  whole  evolution  of  the  linear  perturbations. 

Conclusions 

The  constant-mass  perturbation  model  predicts  that  when  R\  decreases  by  a  factor  of  105, 
ft  increases  by  ~  1010.  This  estimate,  which  is  based  on  simple  geometric  considerations, 
omits  the  dynamical  effects  in  the  problem.  It  is  equivalent  to  setting  the  terms  on  the 
right  hand  side  of  Eq.  (21)  to  zero.  Our  model  predicts  amplification  by  a  factor  of 
approximately  20.  Evidently  the  terms  on  the  right  hand  side  of  Eq.  (21),  though  initially 
small,  play  an  important  role  in  reducing  the  growth  of  the  perturbations.  In  view  of 
the  idealizations  in  our  model  and  of  the  many  assumptions  it  incorporates,  the  coasting 
approximation  is  unlikely  to  be  good  over  a  time  during  which  R\  changes  by  more  than 
10  or  100  and  the  perturbation  amplitude  increases  by  a  factor  of  two  to  three.  Hence  our 
quantitative  results  should  be  regarded  as  being  merely  indicative. 

The  results  of  our  numerical  example  suggest  that  conservation  of  action  is  an  ex¬ 
tremely  reliable  tool  for  predicting  the  evolution  of  surface  waves  on  an  imploding  spherical 
shell.  Although  these  results  are  strictly  applicable  only  to  a  model  in  which  the  shell  is 
uniform,  incompressible,  and  coasting,  the  adiabatic  invariant  (54)  can  probably  be  used 
to  describe  incompressible  irrotational  perturbations  whenever  Ri  <  0  and  the  shell  is 
“thick.”  Note  that  even  though  the  waves  undergo  amplification,  we  arc  not  talking  here 
about  an  instability,  such  as  the  Rayleigh-Taylor  instability.  The  imploding  shell  is  stable 


against  Rayleigh-Taylor  modes;  in  fact,  if  the  sign  of  Ri  were  reversed,  the  present  theory 
would  no  longer  hold. 

The  formalism  which  gave  rise  to  Eqs.  (21)-(24)  can  probably  be  generalized  in  a 
number  of  ways.  First,  the  shell  could  be  made  up  of  several  uniform  layers,  or  even  vary 
as  exp  [/3(i23  —  i2i3)],  where  (3  =  const,  and  still  be  consistent  with  Eq.  (4).  Second,  the 
requirement  of  irrotationality  could  be  relaxed.  Third,  it  is  conceivable  that  additional 
physical  processes  could  be  modeled,  e.g.,  loss  of  material  through  ablation.  The  require¬ 
ment  that  the  basic  state  be  incompressible  is  essential,  however;  otherwise  the  unper¬ 
turbed  motion  of  the  shell  would  be  much  more  complicated  because  the  pressure  buildup 
associated  with  convergence  would  cause  large-amplitude  sound  waves  to  propagate  radi¬ 
ally  through  it.  In  spite  of  this  it  is  plausible  that  incompressible  surface  perturbations  on 
a  compressible  shell  are  accurately  described  by  Eq.  (54).  This  assertion  is  supported  by 
comparison  of  the  behavior  of  the  unstable  case  (the  Rayleigh-Taylor  instability)  in  com¬ 
pressible  and  incompressible  media.  It  has  been  shown2,  e.g.,  that  compressibility  changes 
growth  rates  by  a  less  than  a  factor  of  two. 

In  summary: 

(i)  We  have  defined  the  “coasting  state”  of  an  imploding  spherical  shell  and  derived 
an  analytical  description  of  its  behavior. 

(n)  We  have  derived  a  formalism  for  describing  the  evolution  of  infinitesimal  irrota- 
tional  perturbations  of  an  incompressible  uniform  imploding  spherical  shell. 

(m)  We  have  obtained  the  limiting  form  of  these  equations  for  the  cases  of  a  thin  shell 
and  a  thick  shell,  and  have  exhibited  several  analytical  solutions  for  the  latter. 

(ii>)  We  have  recovered  the  law  of  conservation  of  wave  action  from  the  equation  for 
the  evolution  of  a  surface  wave  on  the  inner  surface  of  a  thick  shell. 

(u)  We  have  numerically  solved  the  exact  equations  for  the  evolution  of  the  coupled 
inner  and  outer  surface  waves  and  shown  that  they  are  accurately  approximated  by  the 
solutions  in  the  thin-  and  thick-shell  limits. 

(m)  We  find  that  for  realistic  compression  factors,  perturbations  of  the  inner  surface 
are  amplified  very  little,  by  only  a  factor  of  order  two  or  three.  This  shows  the  inadequacy 
of  the  simple  picture  in  which  the  mass  of  a  bump  or  asymmetry  per  unit  solid  angle  is 
conserved. 
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RATIO  OF  RADI  I  RAISED  TO  POWER  L 


SHELL  INNER  RADIUS 

Fig.  2  Ratio  (Ri/Ih)1  vs  log  Rx  for  l  =  1,2,5,10,20,50,100,200,500,1000  (reading 
from  right  to  left). 


Fig.  3  Linearly  increasing  amplitude  fi  (or  £2)  of  surface  perturbations  £1 ,  ft  vs  log/?! 
in  the  thin-shcll  approximation,  assuming  initial  unit  amplitudes  and  initially  vanishing 
velocity  coefficients  Ui,it2- 
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AMPLITUDE  FROM  ADIABATIC  INVARIANT 
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SHELL  INNER  RADIUS 

Fig.  4  Amplitude  of  the  inner  surface  oscillations  vs  logi2i,  obtained  from  consci 
tion  of  wave  action.  The  curve  has  been  normalized  to  have  an  initial  value  of  unitv 


PERTURBATION  AMPLITUDE 


INNER  PERTURBATION  VS.  R1  FOR  L  -  1 


Fir.  5(a)  Amplitude  of  inner  surface  wave  vs  log  R\  obtained  from  numerical  solution 
of  the  exact  equations  (2l)-(24)  for  /  =  1  (solid  trace)  and  as  approximated  by  Eq.  (49) 
(broken  trace). 


OUT  L  R  PERTURBATION  VS.  R1  FOR  L  =  1 
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SHELL  INNER  RADIUS 

Fig.  5(6)  Amplitude  of  outer  surface  wave  vs  log  R i  obtained  from  numerical  solution 
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Fig.  6  Amplitudes  of  inner  (solid  trace)  and  outer  (broken  trace)  surface  waves  vs 
log  J?!  for  (a)  /  =  2;  (6)  /  =  5;  (c)  /  =  10;  (d)  l  =  20;  (e)  /  -  50;  (/)  l  =  100;  (g)  l  =  200; 
( h )  l  =  500;  (*)  l  =  1000. 
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Fig.  7  Amplitudes  of  inner  surface  waves  shown  in  Fig.  6,  scaled  by  the  approximate 
function  plotted  in  Fig.  4,  vs  log  .ft  i  for  (a)  l  =  2;  (6)  l  -  5;  (c)  l  =  10;  ( d )  l  =  20;  (e) 
l  =  50;  (/)  l  =  100;  (g)  l  =  200;  (h)  l  =  500;  (»)  l  =  1000. 


